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PREFACE 


This  work  is  related  to  an  ocean-bottom  research  program  at  NRL.  The 
program  is  motivated  by  the  recent  development  of  a  variety  of  bathymetric 
measurement  systems  that  resolve  aspects  of  bottom  variability  previously 
beyond  the  realm  of  quantitative  study.  The  general  aim  of  the  research  is  to 
develop  statistical/stochastic  representations  of  ocean-bottom  variability  at 
scales  that  are  relevant  to  underwater  acoustic  systems.  Such  representations 
will  allow  the  prediction  of  acoustic  system  performance  and  lead  to  the  better 
design  and  use  of  underwater  acoustic  systems  that  either  utilize  or  are  sensi¬ 
tive  to  acoustic  interaction  with  the  ocean  bottom. 


IV 


ANALYTICAL  AND  PHYSICAL  ASPECTS  OF  TWO-DIMENSIONAL 

SPECTRA  ASSOCIATED  WITH 
STATIONARY  RANDOM  PROCESSES 


INTRODUCTION 

This  report  is  an  exploration  cf  certain  fundamental  and  general  aspects  of  spectral  relationships 
that  arise  in  the  analysis  of  random,  or  stochastic,  processes  that  depend  on  more  than  one  dimension. 
The  multidimensional  spectrum  is  the  central  focus  of  our  remarks  rather  than  the  correlation  function. 
A  main  point  of  the  discussion  is  that  several  types  of  spectra  may  arise  in  the  analysis  of  the  multidi¬ 
mensional  situation,  which  can  be  a  source  of  confusion.  There  is  a  significant  difference  with  the 
one-dimensional  situation  which  is  the  case  most  often  treated  in  the  literature.  Stochastic  processes 
that  occur  over  more  than  one  dimension  offer  a  much  greater  richness  in  terms  of  statistical  quantities 
that  can  be  employed  for  their  description.  We  stress  general  aspects  that  may  enter  into  the  considera¬ 
tions  of  many  specific  applications.  Furthermore,  the  physical  interpretation  of  the  various  quantities 
that  occur  is  stressed  as  opposed  to  a  more  formal  mathematical  presentation. 

There  is  a  wide  prevalence  of  power-law  representations  used  to  represent  spectra  in  virtually 
every  field  in  which  spectra  are  used.  This  is  undoubtedly  due  to  the  ease  of  analytically  dealing  with 
such  simple  forms  and  they  further  offer  a  concise  means  of  summarizing  spectral  behavior.  We  evalu¬ 
ate  a  function  that  relates  various  spectral  forms  in  the  case  of  power-law  behavior  and  investigate  its 
analytical  behavior. 

Although  much  of  the  discussion  is  aimed  at  general  properties,  we  give  some  presentation  of  the 
ocean  surface  wave  case  and  suggest  that  it  can  be  viewed  as  a  paradigm  for  the  multidimensional  case. 
This  suggestion  is  partly  motivated  by  the  long  historical  evolution  of  attempts  to  create  a  description  of 
the  dynamics  of  the  sea  surface,  which  has  led  to  a  rich  literature  devoted  to  probabilistic,  or  stochastic, 
models  of  the  sea  surface  [1,2].  Reference  3  discusses  the  different  viewpoints  and  problems  associated 
with  the  early  struggle  to  fashion  a  statistical  approach  to  ocean  wave  spectra.  The  main  point  of  the 
suggestion  is  that  knowledge  of  the  ocean  surface  wave  literature  may  be  a  useful  aid  in  the  description 
and  investigation  of  other  processes. 

VARIOUS  FORMS  OF  SPECTRAL  ENERGY  AND  THEIR  PHYSICAL  INTERPRETATION 

Consider  a  stochastic  process  that  depends  on  two  independent  variables  which  we  designate 
x  andy.  The  primary  question  of  interest  is  the  description  of  the  spectral  structure  associated  with 
7)  (x,y),  where  r\  represents  some  physical  aspect  of  the  process.  Thus  one  can  say  that  evolves  in  a 
random  fashion  as  a  function  of  x  andy.  We  further  assume  that  the  process  is  stationary  [4]  so  that 
the  statistical  structure  does  not  vary  with  location  in  the  plane.  This  assumption  enables  one  to  speak 
in  a  meaningful  way  about  statistical  averages  associated  with  the  process.  Averages  are  denoted  by  an 
overbar.  As  a  matter  of  convenience,  we  assume  that 

V  =  0, 

so  that  we  are  dealing  with  fluctuations  about  a  zero  level. 
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It  is  possible  to  represent  the  process  in  terms  of  a  Fourier-Stieltjes  integral 

t?(x)  =  Jk  dA(k)eikx,  (1) 

where  the  integration  is  over  all  wavenumber  space.  Certain  consequences  of  this  representation  are 
our  main  concern.  The  Fourier-Stieltjes  representation  is  easily  interpreted  as  a  superposition  of  an 
infinite  number  of  infinitesimal,  sinusoidal  waves  of  different  wavenumbers  with  different  directions. 
The  component  associated  with  wavenumber  k  has  an  amplitude  dA  (k)  and  an  orientation  in  physical 
space  specified  by  k.  It  is  convenient  to  describe  some  aspects  of  the  representation  in  terms  of 
wavenumber  space  and  other  aspects  in  terms  of  physical  space  as  indicated  in  Fig.  1. 


(a)  PHYSICAL  PLANE  <b)  WAVENUMBER  SPACE 

Fig.  1  —  Relationship  between  physical  space  and  wavenumber  space  for  the  Fourier-Stieltjes  representation  of  the  sto¬ 
chastic  process.  In  (a)  the  crests  associated  with  the  wave  component  k  are  indicated.  The  crests  are  perpendicular  to 
the  direction  k  and  have  a  spacing  1/1  Ac  I .  The  wave  component  shown  arises  from  a  point  in  wavenumber  space. 


A  fundamental  property  of  the  Fourier-Stieltjes  coefficients  is  that 

dA(k)dA*  (k')  -  0  if  k  5*  k'  (2) 

-  S  (k)  dk  if  k  =  k'. 


From  the  representation  in  Eq.  (1)  it  can  be  shown  that 

^  =  /k  |<//l(k)|2, 

or  using  Eq.  (2), 

=  Jk  S(k)dk.  (3) 

Thus  S( k)  gives  the  density  of  contributions  to  the  mean  square  per  unit  area  of  wavenumber  space. 
S(k)  is  referred  to  as  the  wavenumber  spectrum,  or  power  spectral  density,  of  the  process.  The 
wavenumber  spectrum  tells  us  how  the  mean  square  value  is  distributed  over  the  wavenumber  plane. 


Some  additional  properties  of  the  coefficients  should  be  noted.  Equation  (2)  indicates  that  dif¬ 
ferent  components  are  uncorrelated  as  a  consequence  of  the  stationarity  assumption.  They  are  not 
necessarily  independent.  We  further  assume  that  r)  measures  some  real  valued  aspect  of  the  process, 
and  this  implies  that 


however,  since 


dA  *  (k)  -  dA  (— k); 
S(k)dk  =  \dA  (k)  |2 , 
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we  have 

S(-k)Jk  =  \dA  (— k)  |2  =  |aL4*(k)  |2  =  \dA  (k)  |2  =  5(k) dk , 

so  that 

S(-k)  -  S(k). 


(3) 


One  would  initially  think  that  the  Fourier-Stieltjes  components  should  have  no  constraints  as  far 
as  their  behavior  in  the  wavenumber  plane  is  concerned.  This  is  true  because  such  constraints  produce 
restrictions  on  the  nature  of  the  directionality  of  the  components.  However,  in  the  present  case,  the 
process  17  (x)  consists  of  an  ensemble  of  functions  each  of  which  represents  a  surface  that  is  not  chang¬ 
ing  in  time.  Consequently,  while  the  surface  may  be  represented  in  terms  of  component  sinusoidal 
variations,  the  sinusoidal  components  are  not  propagating,  and  it  is  therefore  impossible  to  distinguish  a 
sinusoid  associated  with  +k  from  one  associated  with  -k.  The  Fourier-Stieltjes  representation  simply 
assigns  an  equal  amplitude  (and  energy)  at  both  wavenumber  k  and  —  k  to  produce  a  real  valued 
sinusoidal  oscillation  of  the  form 


sin  ( kxx  -I-  kyy  +  <j>). 

Later  we  discuss  the  specific  case  of  ocean  surface  wind-waves  where  it  is  possible  to  invoke  a  certain 
assumption  about  the  nature  of  the  process  to  remove  the  apparent  ambiguity. 

This  was  a  summary  of  the  important  properties  of  the  Fourier-Stieltjes  representation  of  the  pro¬ 
cess.  A  mathematical  justification  with  varying  levels  of  rigor  can  be  found  in  several  references  [5-7], 
The  important  quantity  for  the  following  discussion  is  the  spectrum  S(k),  which  can  also  be  written  in 
terms  of  the  wavenumber  components  as 

S(k)  =  S(kx,  ky), 

and  the  relation  of  this  spectrum  to  the  mean  square  tj2  cited  in  Eq.  (3).  We  investigate  several  quanti¬ 
ties  that  are  associated  with  S  (k) .  An  important  ingredient  is  the  physical  interpretation  of  the  quanti¬ 
ties  and  how  they  are  related  to  the  fundamental  spectrum  S( k). 

It  is  sometimes  convenient  to  describe  the  distribution  of  spectral  energy  in  terms  of  polar  coordi¬ 
nates  in  the  spectral  plane.  This  form  would  be  useful  if  one  wanted  to  investigate  the  distribution  of 
spectral  energy  as  a  function  of  angle  in  the  wavenumber  plane  (and  would  correspond  to  directional 
properties  of  the  process  in  physical  space).  The  desired  spectrum  would  have  the  form  S (k,  9 ),  where 
k  is  the  magnitude  of  the  wavenumber  and  9  is  an  angle  in  the  spectral  plane  measured  counter¬ 
clockwise  from  the  kx  axis.  For  this  spectrum  we  have 

V  =  iT/cT  S(k,9)d9dk.  (4) 

We  can  refer  to  S(k,  9)  as  a  directional  spectrum.  It  can  be  thought  of  as  a  distribution  of  mean 
energy  in  k,  9  space  as  indicated  in  Fig.  2  where  a  differential  area  centered  on  k,  9  contributes  an 
amount  S(k,  9)d9dk  to  the  mean  square. 

To  deduce  S(k,  9)  from  S(kx,  ky)  it  is  useful  to  think  of  a  small  differential  angular  sector  in  the 
kx ,  ky  plane  (see  Fig.  2(a)).  Its  contribution  to  the  mean  square  is  S(k,  9)d9dk,  which  can  also  be 
written  as 

S(k,  9)d9dk  =  //  S(kx ,  ky)dkxdky. 

sector 

When  the  differential  sector  is  sufficiently  small  S(kx,  ky)  will  be  effectively  constant  over  the  dif¬ 
ferential  sector  so  that 

S(k,  9)d9dk  =  S(kx,  ky)  Jf  dkxdky  =  S(kx,  ky)kd9dk. 

sector 
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(a) 


(b) 


Fig.  2  —  The  plane  associated  with  S(kxy  ky )  in  (a)  and  that  associ¬ 
ated  with  S(k,  9)  in  (b).  S(ky  0)  arises  from  a  transformation  of  the 
kxy  ky  space  to  k,  9  space.  The  contribution  to  the  mean  square 
from  the  differential  area  in  (a)  is  equal  to  the  contribution  from  the 
differential  area  in  (b). 


This  yields  the  relation 


S(k ,  0)  =  kS(ky  9)  (5) 

in  which  S (k,  9)  designates  the  function  S(kx,  ky)  expressed  in  terms  of  the  polar  coordinates.  The 
relation  between  these  two  spectra  is  thus  established  by  considering  the  transformation  from  one  space 
to  another  or,  equivalently,  from  one  representation  to  another. 


Another  spectrum  of  interest  arises  when  we  describe  the  contribution  to  the  mean  square  made 
by  wavenumber  components  that  are  restricted  to  bands_to  wavenumber  magnitude.  The  situation  is 
indicated  in  Fig.  3.  We  can  define  the  contribution  to  172  from  a  differential  band  of  wavenumbers  as 
S(k)dk  so  that 

r}2  =  J0  S(k)dk.  (6) 

The  spectrum  S  (k)  is  sometimes  referred  to  as  a  scalar  wavenumber  spectrum,  which  is  suggestive  of 
its  origin  as  a  spectral  quantity  that  describes  spectral  energy  irrespective  of  orientation  of  the 
wavenumber.  It  is  often  used  in  the  discussion  of  processes  which  depend  on  more  than  one  spatial 
dimension.  It  is  particularly  useful  in  those  situations  where  the  basic  process  is  isotropic  with  the  fun¬ 
damental  spectrum  S(kx ,  ky)  being  a  function  only  of  the  magnitude  of  the  wavenumber.  For  the  gen¬ 
eral  situation  we  have  that 

S(k)dk  =  //  S(kx ,  ky)dkxdky.  .  (7) 

band 
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Fig.  3  —  Sketch  showing  the  differential  band  in  wavenumber  space  which 
contributes  an  amount  S(k)dk  to  the  mean  square  t/2 

Converting  the  integration  in  Eq.  (7)  to  polar  coordinates  yields 

S(k)dk  =  Jq  '  S(k,  0)kd0dk, 

so  that 

S(k)  =  Jo  S(k ,  0)kde.  (8) 

This  gives  the  relation  of  S(k)  to  the  fundamental  two-dimensional  spectrum  of  the  process. 

The  scalar  wavenumber  spectrum  can  be  thought  of  as  a  one-dimensional  spectrum  in  the  sense 
that  it  describes  the  distribution  of  spectral  energy  over  the  one-dimensional  axis  of  k.  Part  of  its  use¬ 
fulness  stems  from  the  reduction  to  one-dimension  of  the  more  complicated  spectral  energy  distribution 
in  wavenumber  space.  However,  note  that  it  refers  to  a  process  defined  over  two  dimensions. 

It  is  also  possible  to  associate  one-dimensional  spectra  with  the  process  by  asking  for  the  contribu¬ 
tion  to  the  mean  square  made  by  spectral  components  having  a  kx  component  in  a  small  band  of  k 
values.  We  can  write  * 


S(kx)dkx  =  [  f  ^  S(kx,  ky)dkx  ] dky 

so  that  the  one-dimensional  spectral  density  is 

S(kx)  =  J  ^  S(kx,  ky)dky.  (9) 

Figure  4  shows  the  region  of  wavenumber  space  that  contributes  to  the  energy  S(kx)dkx.  The  mean- 
square  is  given  by 

=  /_  S(kx,  ky)dkxdky. 
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Fig.  4  —  Schematic  indication  of  the  region  of  wavenumbers  which  make  a  contri¬ 
bution  to  S(kx)dkx ,  where  S(kx)  is  the  one-dimensional  spectrum  associated  with 
the  x-direction 


which  can  also  be  written  as 

t)2  *=  f  ^  S(kx)dkx  (10) 

in  terms  of  the  one-dimensional  spectrum  S(kx). 


The  symmetry  property  of  S(kx,  ky)  leads  to  a  symmetry  property  of  the  one-dimensional  spec¬ 
trum.  This  can  be  derived  by  evaluating 

S(-kx )  =  S(-kx,  ky)dky. 


From  Eq.  (3) 


S(—kx ,  ky)  =  S(kx,  —ky)  , 


so  that 

S(-kx)  =  J*_“  S(kx,  -ky)dky. 

The  integration  can  be  carried  out  by  changing  variables  to  k'y  =  —ky  to  get 

S(~kx)  =  J_“  S(kx,  k'y)dk'y 

and,  when  compared  to  Eq.  (9) ,  we  conclude 

Si-^-Sik,).  (11) 

This  establishes  the  symmetry  property  for  the  one-dimensional  spectrum. 


An  important  consequence  of  Eq.  (9)  to  remember  is  that  the  one-dimensional  spectrum  S(kx) 
contains  contributions  from  an  infinite  range  of  ky.  This  is  easily  pictured  in  terms  of  Fig.  5,  which 
shows  a  track  parallel  to  the  x-axis  along  which  measurements  of  the  process  are  made  and  from  which 
the  one-dimensional  spectrum  can  be  calculated.  It  is  clear  from  Fig.  5  that  all  spectral  components  of 
the  process  contribute  to  the  variability  of  77 (x,  y )  along  the  sampling  track.  Reference  7  describes  the 
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sampling 

TRACK 


Fig.  5  —  Schematic  illustration  to  demonstrate  how  all  spectral  components  of  the  two- 
dimensional  process  contribute  to  the  variability  along  a  track  parallel  to  the  x-axis.  Each 
of  the  spectral  components  is  defined  throughout  the  plane  and  are  only  partially  indi¬ 
cated  in  the  figure. 


situation  by  terming  the  one-dimensional  spectrum  as  an  "aliased"  version  of  the  two-dimensional  spec¬ 
trum.  Although  we  have  restricted  attention  to  the  one-dimensional  spectrum  S(kx),  it  is  clear  that  an 
infinite  number  of  one-dimensional  spectra  can  be  derived  from  the  process  by  measuring  the  process 
along  lines  at  various  orientations  in  the  physical  plane.  Each  of  the  possible  one-dimensional  spectra 
would  have  the  same  dependence  on  all  components  of  the  process  as  indicated  for  S(kx). 

In  this  section  we  have  considered  a  variety  of  spectra  that  can  be  associated  with  the  random  pro¬ 
cess  r}(x,  y)  and  have  attempted  to  stress  the  physical  interpretation  of  the  spectra  as  well  as  their  rela¬ 
tion  to  the  fundamental  wavenumber  spectrum  S(kx,  ky).  As  the  complexity  of  the  process  increases, 
the  number  of  possible  spectra  increase  and  can  pose  a  serious  problem  to  understanding.  The  problem 
is  compounded  by  the  lack  of  a  standardized  terminology.  It  is  essential  to  keep  in  mind  a  clear 
interpretation  of  the  particular  spectral  quantity  and  its  relation  to  the  complete  process  being  described. 

RELATIONSHIP  OF  THE  ONE-DIMENSIONAL  SPECTRUM  TO  VARIOUS 
TWO-DIMENSIONAL  SPECTRA 

It  is  shown  in  the  previous  section  that  if  the  basic  two-dimensional  spectrum  is  known,  a  variety 
of  other  spectra  may  be  related  to  it.  A  significant  question  to  consider  is  what  one  can  infer  about  the 
two-dimensional  spectrum  from  knowledge  of  the  one-dimensional  spectrum.  For  this  question  there 
appears  to  be  no  answer  that  can  be  stated  in  general  terms  other  than  to  say  that  knowledge  of  the 
one-dimensional  spectrum  places  constraints  on  the  possible  form  of  the  two-dimensional  spectrum. 
We  shall  consider  a  special,  but  important,  case  of  the  question  to  which  a  definite  answer  may  be 
given. 

We  first  start  with  an  expression  giving  the  one-dimensional  spectrum  in  terms  of  the  two- 
dimensional  spectrum  expressed  in  terms  of  polar  coordinates.  We  can  write  Eq.  (9)  as 


(12) 


It  is  convenient  to  change  the  integration  from  one  ky  to  polar  coordinates  using  the  substitutions 


(13) 


when  ky  ^  0  and 


ky  |  A:2  A:*j 


1/2 


(14) 
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when  ky  <  0  (See  Fig.  6.)  Then  we  get 

dky  =  k(k2  —  k2)~x^2dk 

from  Eq.  (13)  and 

dky  =  —  k(k2  —  k2)~xl2dk 

from  Eq.  (14).  These  results  allow  us  to  write  Eq.  (12)  as 

„/f  ,  f°°  kS(k,0)dk  ,  f°°  kS(k,-9)dk  | 

sw ' J*. 

for  kx  ^  0.  For  A:x  <  0,  we  can  use  the  symmetry  property  that  S(—kx )  = 

S(k,  —  6)  in  the  second  integral  arises  because  for  integration  over  the  line  in 
Fig.  6)  we  have  the  relation  that  the  polar  angle  as  a  function  of  k  must  be  the  negative  of  the  polar 
angle  that  appears  in  the  First  integrand. 


(15) 

S(kx).  The  presence  of 
the  fourth  quadrant  (see 


Fig.  6  —  Line  in  the  spectral  plane  along  which  integration  must  be  per¬ 
formed  to  obtain  the  one-dimensional  spectrum  S(kx )  from  the  two- 
dimensional  spectrum  S(kx,  ky).  The  integration  is  most  conveniently 
done  by  considering  the  upper  and  lower  line  segments  separately.  For 
the  upper  segment  9  =  0(k >  kx )  while  on  the  lower  segment  we  can  write 
9  =  —  9(k,  kx).  On  each  line  segment  k  will  vary  from  kx  to  «>. 


It  is  now  useful  to  obtain  a  form  of  the  result  that  applies  when  the  function  S(k,0)  does  not 
depend  on  the  angle  0,  so  that  we  can  write  S(k,9 )  =  S(k).  This  is  the  special  case  when  contours  of 
the  two-dimensional  spectrum  S(kx,ky )  are  circles  in  the  kx,ky  plane  and  arises  in  many  situations  of 
practical  interest.  We  can  refer  to  this  as  the  isotropic  case  and  thereby  define  isotropy  in  terms  of 
spectral  behavior.  A  more  standard  procedure  is  to  define  isotropy  in  terms  of  the  correlation  function 
of  the  process,  but  that  route  would  take  us  too  far  off  the  path  of  the  present  discussion. 
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We  get  from  Eq.  (15)  for  the  case  of  independence  of  angle  the  result  that 


(16) 


where  the  absolute  value  allows  the  expression  to  be  used  for  positive  or  negative  values  of  kx.  The 
quantity  kx  enters  the  integrand  as  a  fixed  parameter  and  is  also  a  limit  on  the  integral.  Equation  (15) 
specifies  the  relationship  between  the  one-dimensional  spectrum  S(kx)  and  the  two-dimensional 
spectrum  S(k ),  which  is  the  two-dimensional  spectrum  in  the  kx,ky  plane  written  as  a  function  of  k. 

We  can  also  establish  a  relationship  between  S(kx)  and  the  isotropic  two-dimensional  spectrum 
S(k).  From  Eq.  (8)  we  have  the  simple  relation  that 


2v  kS(k)  —  S(k). 


(17) 


This  enables  Eq.  (16)  to  be  written  as' 


(18) 


which  gives  an  important  relation  between  the  one-dimensional  spectrum  and  the  isotropic  spectrum 
S(k)  that  refers  to  the  complete  process,  occurring  over  two  dimensions. 

Equation  (18)  gives  rise  to  the  question  as  to  whether  the  relationship  can  be  inverted  so  as  to 
obtain  S(k)  from  a  known  S(kx).  This  problem  arises  in  practice  when  the  process  is  known,  or 
assumed,  to  be  isotropic  but  only  one-dimensional  data  are  available  to  describe  the  process.  The 
inversion  is  made  possible  by  the  theory  associated  with  Abel’s  integral  equation,  of  which  Eq.  (18)  is  a 
particular  example  [8,9].  The  theory  yields  the  inverse  result  that 


(19) 


Thus  for  an  isotropic  process,  the  spectra  S(k )  and  S(kx)  can  be  regarded  as  an  Abel  transform  pair. 
For  some  examples  of  the  formal  use  of  the  Abel  transform  relation  in  various  fields  of  geophysics 
[10,11], 

It  is  desirable  to  give  some  rationale  for  why  there  should  be  such  an  inverse  relationship.  The 
appeal  to  Abel’s  integral  equation  is  essentially  a  formal  mathematical  treatment  of  the  issue  which 
offers  no  intuitive  basis  for  understanding  the  existence  of  the  relationship.  Figure  7  shows  a  useful 
way  to  think  of  the  spectral  plane  for  establishing  an  intuitive  understanding  of  the  inverse  relations. 
The  plane  is  divided  into  circular  bands  of  thickness  dk  which  have  a  mean  energy  content  of  S{k)dk 
and  a  set  of  vertical  bands  of  thickness  dkx  (where  dkx  =  dk)  with  a  mean  energy  content  of  S(kx)dkx. 
It  is  clear  from  the  figure  that  the  circular  bands  establish  a  partition  of  the  vertical  bands.  Further,  the 
mean  energy  of  a  vertical  band  can  be  calculated  by  summing  the  contributions  from  the  various 
regions  of  the  circular  bands  which  make  up  the  partition.  This  will  obviously  involve  geometric  factors 
and  values  of  S(k )  for  k  ^  kx. 

The  following  shows  how  the  inverse  result  of  expressing  S{k)  in  terms  of  S(kx)  seem  reason¬ 
able  (see  Fig.  7).  Consider  the  set  of  values  S(kx )  where  kx  ranges  from  k  to  °°  in  increments  dk. 
Each  value  S(kx)  can  be  represented  as  a  sum  over  S(k)  values,  as  previously  discussed.  Thus  there  is 
a  matrix  relation  between  the  infinite  set  of  S(kx)  values  and  the  infinite  set  of  S(k)  values.  To 
represent  S(k)  in  terms  of  S(kx)  it  is  necessary  to  invert  the  matrix  relationship.  The  result  of  invert¬ 
ing  the  matrix  must  take  the  form  of  Eq.  (19)  when  transition  is  made  from  a  summation  to  an 
integral.  These  remarks  are  intended  to  indicate  a  physical  justification  for  the  inverse  result.  They  do 
not  constitute  a  mathematical  derivation  of  the  result  but  do  suggest  the  lines  along  which  such  a 
derivation  may  be  pursued. 
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Fig.  7  —  Illustration  to  give  rationale  for  expecting  an  inverse 
relationship  to  exist  between  S(kx)  and  S(k).  The  vertical 
bands  contain  an  energy  S{kx)dkx  while  the  circular  bands  con¬ 
tain  an  energy  S(k)dk,  where  dk  =  dkx  for  the  situation  shown. 
The  circular  bands  establish  a  partition  of  the  vertical  bands 
which  makes  it  clear  that  S{kx)dkx  is  the  sum  of  contributions 
arising  from  the  circular  bands.  The  portion  of  each  circular 
band  with  must  be  taken  is  determined  from  geometric  con¬ 
siderations. 


Note  Eqs.  (18)  and  (19)  are  fully  equivalent  to  the  forms  obtained  by  Bell  [10]  starting  with  the 
correlation  function.  To  make  our  equations  conform  with  the  results  of  Ref.  10,  it  is  necessary  to  con¬ 
vert  the  spectrum  S(kx)  to  a  one-sided  form  by  placing  all  the  spectral  energy  into  positive  values  of 
the  wavenumber  kx.  This  would  introduce  a  factor  of  two  and  make  our  results  identical  with  Bell’s. 

It  is  useful  to  summarize  the  relations  that  are  known  to  exist  between  the  different  spectral 
quantities.  This  is  done  in  Table  1  for  the  general  case  and  in  Table  2  for  the  special  case  of  isotropy. 
An  important  aspect  of  Table  1  is  that  some  of  the  relations  are  not  known  in  general  because  of  the 
nonexistence  of  appropriate  inverse  results.  Table  2,  on  the  other  hand,  has  explicit  formulas  for  all 
the  entries,  this  is  entirely  a  consequence  of  the  assumption  of  isotropy.  The  situation  becomes  more 
involved  with  processes  that  depend  on  more  than  two  independent  variables,  since  spectra  could  be 
defined  to  account  for  several  different  directional  aspects  of  the  spectral  energy.  In  such  cases  it  is 
helpful  to  keep  in  mind  the  matrix  of  possible  relations  corresponding  to  Tables  1  and  2  for  the  two- 
dimensional  case. 
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Table  1  —  Relations  between  the  spectral  quantities  for  the  general  case.  The  table  should  be 
read  as  the  quantity  appearing  at  the  left  of  a  row  and  column  entry  is  equal  to  the  specified 
function  of  the  quantity  that  appears  at  the  top  of  the  column.  Note  that  a  column  is  the  set  of 
inverse  relations  for  entries  appearing  in  the  row  having  the  same  label  as  the  column.  Some 
entries  indicate  no  general  relation  because  an  inverse  relation  is  needed  for  which  no  general 
result  is  known. 


S(kx,ky) 

S(k,0) 

i(k.0) 

i(k) 

S(kx) 

S(kx,ky) 

S(k,0) 

= 

CD 

1C 

1-I-* 

NO  GENERAL 

RELATION 

NO  GENERAL 

RELATION 

,$(k,0) 

kS(k,0) 

= 

NO  GENERAL 

RELATION 

NO  GENERAL 

RELATION 

i(k) 

2t r 

1  S(k,Q  )kd0 

0 

2  7T 

/i(k,0)d0 

0 

= 

NO  GENERAL 

RELATION 

S(kx) 

1 S  (kx,ky)dky 

7i(k,0)dk 

4  h‘-k- 

+  7  ^  (k,-  0)dk 
K\  'I  k2-kx 

NO  GENERAL 

RELATION 

= 
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Table  2  -  Relations  between  the  spectral  quantities  for  the  special  case  of  isotropy,  for 
which  S(kx,ky)  may  be  written  entirely  as  a  function  of  k.  This  function  is  designated 
as  S(k),  which  represents  a  change  in  notation  from  Table  1.  The  other  quantities  are 
designated  as  in  Table  1.  The  table  should  be  read  as  the  quantity  appearing  at  the  left 
of  a  row  and  column  entry  is  equal  to  the  specified  function  of  the  quantity  that  appears 
at  the  top  of  the  column  Note  that  a  column  is  the  set  of  inverse  relations  for  entries 
appearing  in  the  row  having  the  same  label  as  the  column  and  that  all  the  inverse  rela¬ 
tions  are  known  for  the  isotropic  case. 


S(kXlky) 

S(k) 

^(k,e> 

^<k> 

S(kx) 

S(kx,ky) 

S(k) 

= 

£i<k,0) 

1  j(k) 

2TTk 

1  7  1  dS(kx) 

nJ  'Tkl^k7’  dkx  dK 

k 

^(k,0) 

kS  (k) 

= 

i^<k) 

1  7  k  ds(kx)  j( 
"J'JkZ-k2  dkx 

k 

i(k) 

2  tt  k  S(k) 

2ttJ)  (k,0) 

— 

7  k  dS(kx)  „ 

J  4 kx-k2  dkx 

S(kx) 

oo 

[ S(kx,ky)dky 

^7  ^(k,0)dk 

A  'Jk!-k- 

1  7  ^  (k)dk 

fk2-k2 

M  VK  k* 

EVALUATION  OF  THE  SPECTRA  FOR  POWER-LAW  MODELS 

An  important  and  useful  analytical  representation  of  spectra  is  offered  by  power-law  approxima¬ 
tions  to  empirical  spectra.  This  is  especially  true  for  physical  processes  that  are  of  a  wide-band  nature 
and  for  which  the  spectral  energy  is  distributed  over  the  band  with  no  strong  preference  for  isolated 
wavenumbers. 

In  particular,  for  the  two-dimensional  wavenumber  spectrum  S(kx,  ky)  it  is  sometimes  useful  to 
use  the  representation 


S(kx,  ky)  —  S0k  “,  (20) 

where  one  might  in  general  assume  that  Sq  and  a  are  functions  of  the  angle  6  in  the  spectral  plane. 
This  type  of  representation  would  be  useful  in  those  situations  where  a  strict  power-law  behavior  may 
be  applicable  along  individual  rays  from  the  origin  but  with  the  parameters  of  the  power  law  depending 
on  which  particular  ray  is  being  used. 
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Little  additional  simplification  can  be  achieved  with  the  expression  in  Eq.  (20)  when  S0  and  a  are 
allowed  to  vary  with  9.  Now  consider  the  special  case  of  isotropy  for  which  the  two-dimensional  spec¬ 
trum  S(kx,  ky)  can  be  represented  as  in  Eq.  (20)  with  constant  values  of  S0  and  a. 

The  scalar  wavenumber  spectrum  is  given  by  Eq.  (8)  as 

S(k)  -  2vkS , 


or 

S(k )  =  2tt  50A:-a+1.  (21) 

An  important  aspect  of  the  scalar  wavenumber  spectrum  S(k)  is  that  the  exponent  in  the  power  law  is 
different  from  the  exponent  which  appears  in  the  expression  for  S(kx,  ky ). 

Another  spectrum  of  interest  is  the  one-dimensional  spectrum  S(kx )  that  would  arise  in  a  study 
of  behavior  of  the  process  tj(x,  y)  along  the  x-axis.  Since  the  process  is  isotropic,  there  is  no  special 
significance  attached  to  the  x-axis;  the  statistical  behavior  of  the  process  along  any  line  in  space  will  be 
identical  with  the  behavior  along  any  other  line  in  space.  There  are  several  ways  to  obtain  S(kx)y  but  it 
is  perhaps  simplest  to  derive  it  from  the  basic  relation 

S{kx)  =  S(kx,  ky)dky. 

Using  the  form  in  Eq.  (20)  and  the  symmetry  in  ky,  this  becomes 

S(kx )  =  2  Jo°°  S0k~adky.  (22) 

The  evaluation  of  this  integral  is  carried  out  by  considering  kx  as  positive  and  using  the  substitutions 

k  »  y/itx  +  ky* 
ky  =*  kx  tan  9 , 
dky  —  kx  sec 2  9d9 . 

The  result  is  that  Eq.  (22)  becomes 

3L 

S(kx)  =  2S’0|Arx|“o,+l  J^2  cos a~20de.  (23) 

The  one-dimensional  spectrum  is  written  in  a  form  applicable  to  positive  or  negative  kx .  It  has  a 
power-law  behavior  that  differs  from  that  of  the  two-dimensional  spectrum  by  having  an  exponent  that 
is  increased  by  1.  The  situation  is  quite  similar  to  the  relationship  between  the  two-dimensional  spec¬ 
trum  and  the  scalar  wavenumber  spectrum  S (k).  The  only  difference  between  the  mathematical 
expressions  for  S(kx )  and  S(k)  involves  a  constant  factor  that  arises  because  these  two  spectra  refer  to 
different  types  of  energy  associated  with  the  process.  This  difference  is  indicated  in  the  section  on  vari¬ 
ous  forms  of  spectral  energy. 

It  is  desirable  to  evaluate  the  one-dimensional  spectrum  explicitly  for  any  situation  that  might 
arise  in  applications.  To  do  this  we  write  Eq.  (23)  as 

S(kx)  =  2 S0\kx\-a+1f(<*  -  2),  (24) 

where 

r»7r/2 

f(a)  =  I  cos  a9dG.  (25) 

*/0 

Now  consider  properties  of  the  function  /(a)  and  the  evaluation  of  the  function  for  specific  values  of 
a. 


The  first  question  to  ask  is  the  range  of  values  for  which  the  integral  in  Eq.  (25)  makes  sense. 
For  values  a  <  —  1  the  integrand  in  Eq.  (25)  becomes  infinite  at  -y  in  such  a  way  that  /(a)  would  be 
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ANGLE  0,  radians 

Fig.  8  —  Behavior  of  the  function  cosa0  for  various  values  of  a 


infinite.  Consequently,  it  is  only  meaningful  to  consider  values  of  a  >  —  1.  Figure  8  shows  an  indica¬ 
tion  of  the  behavior  of  the  integrand  in  Eq.  (25). 

One  important  general  property  of  the  function  f(a)  is  that  it  is  a  decreasing  function  of  a.  This 
is  seen  by  evaluating  the  derivative  of  the  function.  We  get 

d  dr 77  ^ 

——/(a)  =  ——  I  cos  a9d6 

da  da 

(•■nil 

=  Jo  cos“0  log  cos  9  d9.  ,  (26) 

The  integrand  in  Eq.  (26)  is  negative  except  possibly  at  an  endpoint  of  the  integration  interval.  Conse¬ 
quently,  we  have 

~T~f(a)  <  0, 
da 

which  establishes  the  decreasing  nature  of  the  function  in  the  interval  a  >  —  1. 

For  integer  values,  the  function  /(a)  can  be  evaluated  directly  by  a  repeated  application  of  the 
technique  of  integration  by  parts.  This  yields  the  result  that 

r 17  ^ 

f(2n)  =  cos2w  OdO 


2  n  J0 

=  —  1  2n  —  3  1_  v ■_ 

In  2n  —  2  2  2 ’ 
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where  there  are  n  factors  (m  ^  0)  preceding  the  last  factor  of  7r/2.  The  result  for  an  odd  value  of  a  is 
derived  similarly  and  we  get 

/•tt/2 

/(2m  +  1)  =  I  cos2n+  Gd0 
*/0 


In 

2n  +  1 


p  tt!  2 

Jo 


cos2”  l9d9 


In  In  —  2  2_  (28) 

2/2  +  12/1  —  1  3 


where  w  factors  («  >  0)  appear  on  the  right-hand  side  of  Eq.  (28). 


Some  general  relations  can  be  established  between  /  (a),  the  beta  function,  and  the  gamma  func¬ 
tion  which  can  be  used  for  evaluating  f(a)  at  noninteger  values  of  the  argument.  An  integral 
representation  of  the  beta  function  is 


B(x,y)  =  2  f  2  sin2jc  1  9  cos2y  1  9  d9 

where  x,y  are  positive.  The  beta  function  is  related  to  the  gamma  function  by 

T(x)  T(y) 


where  the  gamma  function  satisfies 


r(*  +  ,) 
rOc  +  l)  =  x  r  (x). 


(29) 

(30) 

(31) 


A  comparison  of  the  integral  expression  for  /  (a)  with  the  integral  representation  for  the  beta  function 
in  Eq.  (29)  leads  to  the  result  that 


-L  a  *b  1 
2’  2 


(32) 


which  can  be  expressed  in  terms  of  the  gamma  function  as 

r 

/(.)  -  ds.  - 

r 


a  +  1 
2 


a  +  2 
2 


(33) 


The  notation  and  properties  of  the  beta  and  gamma  functions  are  given  in  Ref.  12. 


Equation  (33)  provides  the  most  useful  means  for  evaluating  /(a)  for  general  values  of  a.  One 
can  simply  use  Eq.  (31)  in  conjunction  with  tabulated  values  of  the  gamma  function  for  some  standard¬ 
ized  interval.  Figure  9  is  a  plot  of  the  function  f(a)  over  a  limited  range  of  a.  It  agrees  with  properties 
already  derived,  from  the  above  results.  The  function  f(a)  is:  (a)  positive-valued,  (b)  tends  to  infinity 
at  a  =  —  1,  (c)  steadily  decreases  in  value  as  a  increases,  (d)  changes  slowly  in  value  when  a  >  >1, 
and  (e)  asymptotically  approaches  zero  as  a  approaches  infinity. 


COMMENTS  ON  OCEAN  SURFACE  WIND-WAVE  SPECTRA 

An  interesting  and  instructive  example  of  the  use  of  spectra  for  a  process  that  depends  on  several 
variables  arises  in  the  study  of  ocean  surface  waves  [1,2].  A  strong  motivation  for  consideration  of  this 
case  comes  from  the  fact  that  ocean  surface  waves  are  a  common  aspect  of  our  experience  with  at  least 
some  aspects  of  the  process  readily  observable.  We  shall  use  ocean  surface  waves  as  a  specific  geo¬ 
physical  model  to  illustrate  some  important  general  points  concerning  spectral  analysis  for  random 
processes  dependent  on  several  variables. 
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Fig.  9  —  The  behavior  of  the  function  f(a)  which  is  related  to  the 
one-dimensional  spectrum  through  S(kx )  =  2S0\kx  |-a+1/(a  -  2) 


The  theoretical  investigation  of  the  spectral  characteristics  of  ocean  surface  waves  uses  a  represen¬ 
tation  of  the  elevation  of  the  ocean  surface  above  some  mean  level  given  by  a  Fourier-Stieltjes  integral, 


dA  (k ,m)  exp 


/  (k  *  x  —  mt) 


(34) 


where  the  integration  is  over  all  wavenumber,  frequency  space.  This  representation  is  obtained  on  the 
basis  that  the  process  may  be  considered  as  a  stationary  random  function  of  both  position  and  time.  An 
important  consequence  of  the  stationarity  is  that  the  Fourier-Stieltjes  coefficients  satisfy 


dA  (k ,m)dA  (k',m')  —  0  if  k,m  ^  k,m' 

=  X (k,m) dkdm  if  k  =  k\m  —  m.  (35) 

It  can  be  shown  that 

7}2  =  f  f  X(k,m)dkdm  (36) 

so  that  X(k,m)  is  the  wavenumber-frequency  spectrum. 

If  it  is  assumed  that  the  ocean  wave  system  consists  of  individual  wave  components  whose  ampli¬ 
tudes  are  sufficiently  small,  then  in  deep  water  we  have  the  dispersion  relation 

O-2  =  gk  (37) 

between  the  frequency  cr  of  an  individual  component  and  the  magnitude  of  its  wavenumber.  A  conse¬ 
quence  of  the  dispersion  relation  is  that  spectral  energy  in  wavenumber-frequency  space  can  only  exist 
in  certain  regions.  Specifically,  we  must  have 

X(k ,  m)  —  i/>(k)8(/w  —  cr),  (38) 
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where  a  is  the  function  of  k  given  in  Eq.  (37).  When  this  result  is  substituted  in  Eq.  (36),  the  result 
is  that 

~2  =  (k)rfk,  (39) 

where  t/»(k)  is  the  wavenumber  spectrum. 

Our  main  reason  for  considering  the  ocean  wave  case  is  that  it  demonstrates  some  interesting 
aspects  of  spectral  analysis.  For  example,  suppose  that  one  made  detailed  measurements  of  the  ocean 
surface  elevation  at  one  instant  of  time  over  a  sufficiently  large  area  of  the  surface.  Such  data  would 
allow  an  estimate  of  the  wavenumber  spectrum  (k),  where  the  result  would  be  obtained  that 

t  (k)  =  ■/,  (-k).  (40) 

That  is,  the  spectral  energy  appearing  at  k  is  the  same  as  at  — k  (see  Fig.  10).  An  actual  example  of  the 
wavenumber  spectrum  derived  from  field  data  is  shown  in  Ref.  1.  The  reason  for  the  equal  distribution 
of  spectral  energy  at  k  and  — k  is  due  to  the  limited  information  assumed  available  for  the  process.  An 
instantaneous  snapshot  of  the  sea  surface  does  not  allow  one  to  distinguish  between  waves  traveling  in 
the  +k  and  -k  directions.  Additional  information  about  the  process  would  be  needed  to  fully  describe 
the  propagational  characteristics  of  the  wave  energy.  This  exemplifies  the  important  point  that  the  par¬ 
ticular  spectrum,  or  spectral  characteristics,  one  can  study  depends  on  the  nature  of  the  measurements. 


Fig.  10  —  Schematic  illustration  of  the  anisotropic  nature  of  the 
wavenumber  spectrum  S(k )  for  ocean  surface  wind-waves.  The  indicated 
spectrum  is  representative  of  what  would  be  calculated  from  a  "snapshot" 
of  the  ocean  surface.  The  spectrum  has  a  symmetry  S(—  k)  =  S(k)  as  a 
result  of  the  type  of  data  assumed,  so  any  half-plane  will  be  an  image  of 
the  other  half-plane. 

Perhaps  the  most  significant  aspect  of  the  spectrum  is  its  obvious  lack  of  isotropy.  The  variation 
of  the  spectral  energy  density  along  a  ray  from  the  origin  will  be  found  to  change  with  direction,  or  the 
particular  ray  being  examined.  Another  aspect  of  such  spectra  is  that  the  spectrum  does  not  continually 
increase  as  the  wavenumber  is  reduced.  This  particular  pattern  of  behavior  is  unusual  for  geophysical 
processes,  which  often  demonstrate  the  characteristics  of  a  red-noise  process. 
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Theoretical  spectra  used  for  the  description  of  ocean  surface  waves  are  often  constrained  in  such  a 
way  as  to  require  that  spectral  components  be  traveling  in  the  downwind  direction.  This  particular  con¬ 
straint  arises  from  observational  studies  of  surface  waves  as  well  as  from  energetic  considerations  that 
identify  the  wind  as  the  primary  source  of  energy  for  the  surface  waves.  The  restriction  can  be  put  in 
the  form  that  nonvanishing  spectral  energy  should  only  appear  at  wavenumbers  with 

k  •  U  ^  0, 

where  U  is  the  wind  vector  at  some  canonical  height  above  the  ocean  surface.  This  allows  the 
wavenumber  spectrum  to  be  represented  as  in  Fig.  11,  where  there  is  now  no  ambiguity  in  the  spectral 
plane.  Theoretical  surface  wave  spectra  are  usually  represented  so  that  a  symmetry  in  the  spectrum 
occurs  relative  to  the  wind  direction,  which  is  also  indicated  in  the  figure.  An  immediate  implication  of 
the  figure  is  that  the  natural  coordinate  system  to  employ  in  studies  of  the  two-dimensional  statistical 
structure  of  ocean  surface  waves  would  have  one  axis  aligned  with  the  wind  direction.  The  existence  of 
a  preferred  direction  oriented  with  the  wind  direction  also  occurs  in  more  general  spectral  wave  models 
[131. 


Fig.  1 1  —  Schematic  illustration  of  the  wind-wave  spectrum  under  the 
assumption  of  down-wind  propagating  waves 


The  importance  of  these  aspects  of  ocean  surface  waves  is  that  they  exemplify  certain  features  of 
random  processes  that  vary  over  a  multidimensional  space  that  must  be  dealt  with  in  other  situations 
not  as  well  understood.  Consider  the  possible  existence  of  an  orientation  direction  for  such  processes. 
Orientation  does  not  arise  as  an  issue  in  truly  one-dimensional  processes.  We  here  argue  that  aniso¬ 
tropy  must  be  dealt  with  on  an  equal  footing  with  stationarity  in  the  treatment  of  data.  The  reason  for 
this  is  that  in  the  collection  of  several  separate  samples  to  form  an  ensemble,  the  data  samples  must  be 
oriented  in  the  same  way  relative  to  one  another  for  ensemble  averages  to  properly  represent  statistical 
structure  associated  with  the  process. 

An  indication  of  the  problem  can  be  obtained  by  thinking  of  the  surface  wave  case.  Suppose  one 
collected  finite  two-dimensional  samples  of  the  ocean  surface  wave  elevation  from  many  areas  of  the 
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world’s  oceans.  These  samples  could  be  used  to  construct  an  ensemble  as  indicated  in  Fig.  12.  Sup¬ 
pose  that  only  samples  with  the  same  fetch,  duration,  and  wind  speed  were  placed  in  the  ensemble  with 
the  additional  restriction  that  all  samples  would  pertain  to  deep  water.  We  are  artificially  creating  an 
ensemble  in  which  only  the  wind  direction  has  not  been  properly  considered.  The  ensemble  average 
spectrum  for  the  ensemble  indicated  in  the  figure  would  turn  out  to  be  in  error.  In  effect,  the  correct 
spectrum  indicated  in  Fig.  1 1  would  be  spread  over  the  spectral  plane  and  much  reduced  in  amplitude. 
This  effect  would  entirely  be  the  result  of  the  failure  to  align  the  data  samples  with  the  same  orienta¬ 


tion. 


Fig.  12  —  Schematic  illustration  of  the  nature  of  ensemble  in  which 
direction  was  not  taken  account  of  properly 


Although  the  above  example  is  rather  artificially  contrived,  it  does  dramatize  the  serious  problems 
associated  with  orientation.  For  a  process  dependent  on  several  variables  and  for  which  little  is  known 
concerning  orientation,  it  is  clear  that  techniques  are  needed  for  determining  orientation  properties  of 
the  process.  There  is  also  a  practical  necessity  for  being  able  to  determine  orientation  in  the  earliest 
stages  of  the  data  processing  if  an  ensemble  type  approach  is  to  be  used.  This  problem  does  not  seem 
to  have  been  adequately  addressed  in  the  literature  on  stochastic  processes. 

There  are,  then,  certain  general  aspects  of  ocean-surface-wave  spectra  that  are  useful  in  a  concep¬ 
tual  way  for  the  study  of  other  processes.  The  main  features  of  the  ocean-surface  wavenumber  spec¬ 
trum  are  the  existence  of  a  preferred  direction  and  the  location  of  the  maximum  in  the  spectral  energy 
density  at  a  wavenumber  whose  magnitude  changes  with  wind  conditions.  The  description  of  ocean 
surface  waves  as  a  stochastic  process  is  in  a  fairly  advanced  state  relative  to  the  description  of  other 
geophysical  processes  and  represents  an  example  of  a  process  dependent  on  several  dimensions  which  is 
useful  for  dramatizing  certain  general  issues  associated  with  such  processes. 

CONCLUSIONS 

We  have  considered  a  variety  of  spectra  which  arise  in  the  discussion  of  multidimensional  stochas¬ 
tic  processes  with  particular  emphasis  on  the  two-dimensional  situation.  The  number  of  spectra  that 
can  be  used  in  the  description  and  analysis  of  such  processes  is  large  and  can  be  a  source  of  confusion. 
We  have  stressed  the  need  to  keep  in  mind  the  physical  interpretation  of  the  different  spectra,  which 
are  related  to  different  ways  of  describing  the  spectral  energy.  The  basic  reason  for  the  large  number  of 
possible  spectral  quantities  for  such  multidimensional  processes  is  that  it  is  possible  to  view  the  process 
in  different  ways.  A  point  worth  stressing  is  the  necessity  to  keep  the  full  nature  of  the  process  in 
mind. 
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The  various  spectra  that  are  introduced  are  all  related  to  the  basic  spectrum  S(kx,  ky)  of  the  pro¬ 
cess  7](x,  y).  A  particular  question  of  interest  is  establishing  inverse  relationships  where  possible.  An 
inverse  relationship  of  particular  importance  is  that  between  the  one-dimensional  spectrum  and  the 
two-dimensional  isotropic  spectrum  in  the  case  of  an  isotropic  process  for  which  the  two-dimensional 
spectrum  is  independent  of  angle  in  the  spectral  plane.  This  is  done  through  appeal  to  Abel’s  integral 
equation,  which  indicates  that  S(kx)  and  S(k)  may  be  thought  of  as  an  Abel  transform  pair.  A  physi¬ 
cal  justification  is  given  for  the  existence  of  the  inverse  relation  to  provide  an  intuitive  understanding 
of  how  the  Abel  transform  relation  arises.  The  physical  justification  is  in  terms  of  an  examination  of 
two  ways  of  describing  the  distribution  of  spectral  energy  in  the  spectral  plane.  It  represents  a  some¬ 
what  novel  way  of  thinking  of  the  inverse  relation  and  seems  capable  of  being  put  on  a  firm  mathemati¬ 
cal  basis. 

We  consider  the  important  case  of  power-law  models  for  the  two-dimensional  spectrum  because  of 
its  expected  use  in  representing  two-dimensional  spectral  behavior.  The  power-law  representation  is 
mainly  considered  for  isotropic  processes.  Again,  the  one-dimensional  spectrum  requires  some 
involved  calculations.  The  one-dimensional  spectrum  is  related  to  a  definite  integral  containing  the 
exponent  of  the  power-law  as  a  parameter.  Some  general  properties  are  established  for  the  integral.  An 
explicit  evaluation  is  given  in  terms  of  the  gamma  function.  The  definite  integral  is  also  plotted  as  a 
function  of  the  parameter  to  provide  a  useful  means  for  obtaining  values  in  particular  cases  which 
might  arise  in  applications. 

The  representation  of  ocean  surface  wind-waves  as  a  stationary  process  is  considered  because  the 
wind-wave  case  offers  a  useful  example  of  a  process  taking  place  over  a  multidimensional  space.  Its 
utility  stems  from  the  rather  extensive  work  done  on  the  problem  and  because  of  the  ease  of  observing 
some  aspects  of  the  process,  even  during  routine  daily  life.  The  windwave  case  is  also  important 
because  it  serves  as  an  example  of  strong  anisotropy  of  the  two-dimensional  spectrum  relative  to  the 
direction  of  the  wind.  A  general  implication  of  the  discussion  of  wind-waves  is  that  anisotropy  of  the 
process  must  be  considered  equally  important  with  stationarity  in  the  treatment  of  data,  especially  if 
ensemble  averages  are  to  be  used  to  estimate  the  two-dimensional  spectrum.  This  is  due  to  the  need  to 
properly  orient  the  different  realizations  that  make  up  the  ensemble  prior  to  averaging.  A  failure  to 
properly  account  for  orientation  can  lead  to  ensemble  average  spectra  that  grossly  misrepresent  the  sta¬ 
tistical  structure  of  the  process. 
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GLOSSARY 


dA  (k ),dA  (k,m) 

Fourier-Stieltjes  coefficients 

B(x,y) 

beta  function 

f(oc) 

special  function  entering  spectral  calculation 

g 

acceleration  of  gravity 

k 

vector  wavenumber 

kx,ky 

components  of  wavenumber 

m 

frequency 

50 

parameter  entering  power  law 

S(k),S(kx,ky) 

two-dimensional  wavenumber  spectrum 

S(k,0) 

two-dimensional  wavenumber  spectrum  expressed  in  polar  coordinates 

s(k,e ) 

directional  spectrum 

S(k) 

scalar  wavenumber  spectrum 

S(kx) 

one-dimensional  spectrum 

u 

wind  vector 

X 

position  vector 

x,y 

spatial  coordinates 

X(k,m) 

wavenumber-frequency  spectrum  for  wind-waves 

a 

exponent  of  power-law 

T(x) 

gamma  function 

80c) 

Dirac  delta  function 

V 

the  fundamental  variable  of  the  stochastic  process  of  interest,  taken  to  be  real- valued 

Zd- 

V 

mean  square  of  19 ,  sometimes  more  casually  referred  to  as  the  mean  energy 

e 

angle  in  spectral  plane 

C T 

frequency 

4> 

phase 

i//(k) 

wavenumber  spectrum  for  wind-waves 

0 

average  of  ( ) 

()* 

complex  conjugate  of  ( ) 
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